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1 Introduction 

In this paper, we investigate the properties of traveling wave solutions to hyper- 
bolic conservation laws augmented with diffusion and dispersion, and reviews 
the existence and qualitative properties of the associated kinetic functions. Such 
a function characterizes the class of admissible shock waves, both compressive 
and undercompressive, selected by a given regularization. Building on the pi- 
oneering papers [24l,l25'TrJi6], the mathematical research on undercompressive 
shocks generated by diffusive-dispersive limits developed intensively in the 
last fifteen years. For backgrovind on this topics and further material, we refer 
the reader to the reviews [|T7l[T8l[T9il and the extensive literature cited therein. 
The present review restrict attention to traveling waves and to a class of scalar 
equations. 

The kinetic relation can be defined as follows. Recall that classical com- 
pressive shocks with a given left-hand state U- (and wave family, when sys- 
tems of equations are considered) form a one-parameter family of solutions, 
parametrized by their right-hand state m+. By contrast, given any left-hand 
state U- (and wave family), there typically exists a single undercompressive 
shock, and the kinetic function (p^ precisely determines the right-hand state 

u+ = (p\u-) 

as a function of the left-hand side. 

The fundamental questions of interest are the following ones: do there 
exist traveling wave solutions associated with classical and/or with nonclassical 
shock waves ? Can one associate a kinetic function to the given model ? If so, 
is this kinetic function monotone ? What is the behavior of arbitrarily small 
shocks ? How does the kinetic function depend upon the parameters? 

Answers to these questions were obtained first for the cubic flux function, 
by deriving explicit formulas for the kinetic function in Shearer et al. [15J and 
Hayes and LeFloch [12J. General flux-functions and general regularization 
were covered by Bedjaoui and LeFloch in the series of papers [3J-[7]. 

More generally, the existence and properties of traveling waves for the 
nonlinear elasticity and the Euler equations are known in both the hyperbolic 
[22, 4 1 and the hyperbolic-elliptic regimes [25 23. 8, 5|. For all other models, 
only partial results on traveling waves are available. 

The existence of nonclassical traveling wave solutions for the thin liquid 
film model is proven by Bertozzi and Shearer in |10J. For this model, no 
qualitative information on the properties of these traveling waves is known, 
and, in particular, the existence of the kinetic relation has not been rigorously 
established yet. The kinetic function was recently determined numerically in 
LeFloch and Mohamadian [20J. For the 3x3 Euler equations, we refer to |[7| . 

Finally, we also recall that the Van de Waals model admits two inflection 
points and leads to multiple traveling wave solutions. Although the physical 
significance of the "second" inflection point is questionable, given that this 
model is extensively used in the applications it is important to investigate 



whether additional features arise. Indeed, it is established in (21 that non- 
monotone nonclassical traveling wave profiles exist, and that a single kinetic 
function is not sufficient to single out the physically relevant solutions. 

An outline of this paper is as follows. In Section 2, we briefly discuss 
the case of the diffusion model, while the rest of the paper is concerned with 
the diffusion-dispersion model. We then begin with the case of a cubic flux- 
function for which explicit formulas can be derived. The main results are stated 
in Section 4 for general flux-functions having one inflection point. Sections 5 
and 6 are concerned with the derivations of key properties of the traveling 
waves and kinetic function, corresponding to a fixed shock speed and to a fixed 
diffusion over dispersion ratio, respectively. 

2 Traveling waves associated with the nonlinear dif- 
fusion model 

Consider the scalar conservation law 

dtu + d:,f{u) = Q, u = u{x, t) e ]R, (2.1) 

where / : ]R ^ IR is a smooth mapping. We begin, in this section, with the 
nonlinear diffusion model 

dtu + (9.v/(w) = £ {b{u) M;c) , u = u'ix, t) e R, (2.2) 

where £ > is a small parameter. The diffusion function b :^ ^ 1R+ is assumed 
to be smooth and bounded below: 

b(u) > fo > 0, (2.3) 

so that the equation (|2.2|| is uniformly parabolic. We are going to establish that 
the shock set associated with the traveling wave solutions of (|2.2b coincides 
with the one described by Oleinik entropy inequalities (see l|2.9t , below). 

Recall that a traveling wave of l|2.2|l is a solution depending only upon the 

variable 

X- At ,^ ^, 

y - ^— (2.4) 

for some constant speed A. Note that, after rescaling, the corresponding trajec- 
tory y 1-^ u{y) is independent of the parameter e. Fixing the left-hand state U- 
we search for traveling waves of l|2.2b connecting^ m_ to some state u+, that is, 
solutions y i-^ u{y) of the ordinary differential equation 

- A My + f{u)y = (b{u) Uy) (2.5) 

satisfying the boundary conditions 

lim u{y) = U-, lim u{y) = u+, lim Uy{y) = 0. (2.6) 

y^-co y— > + oo |y|_j + co ■ 



In view of II2.6I I the equation jl.Si can be integrated once: 

b{u{y))Uy{y) = -Aiuiy)-iQ+fiu{y))-f{u.), yeU. (2.7) 

The Rankine-Hugoniot condition 

-A(u+-m_)+/(m+)-/(m_) = (2.8) 

follows by letting y —> +00 in ||2.7| |. The equation i2.7i is an ordinary differential 
equation (O.D.E) on the real line. The qualitative behavior of the solutions is 
easily determined, as follows. 

Theorem 2.1 (Diffusive traveling waves). Consider the scalar conservation law 
l|2.1|l with general flux-function f together with the diffusive model l|2.2) l. Fix a left- 
hand state U- and a right-hand state u+ + m_. Then, there exists a traveling wave of 
l|2.7|l associated with the nonlinear diffusion model (|2.2b if and only ifu- and u+ satisfy 
Oleirdk entropy inequalities in the strict sense, that is: 

f{v)-f{u-) f{u+)-fiu-) 

> for all V lying strictly between U- and u+. (2.9) 

v-u- u+ - u- -^ •' ^ -^ 

Proof. All the trajectories of interest are bounded, i.e., cannot escape to infinity. 
Namely, the shock profile satisfies the equation 

^^,_ u-u. I f{u)-f{u.) f{u,)-f{u.) \ ^2 

b{u) \ u - U- u+ - u- I ' 

It is not difficult to see that the solution exists and connects monotonically U- 
to M+ provided Oleinik entropy inequalities hold and the right-hand side of 
l|2.10b keeps (strictly) a constant sign (except at the end point y = ±00 where it 
vanishes). D 

We define the shock set associated with the nonlinear diffusion model as 
S(m-) := |m+ / there exists a solution of (|2.6l l — l|2.81 l|. 

From Theorem l2.1l one can deduce the following. 

Theorem 2.2 (Shock set based on diffusive limits). Consider the scalar conservation 
law (|2.1) l when the flux f is convex, concave-convex, or convex-concave (see l|4.2) l, 
below). Then, for any m_, the shock set S{u-) associated with the nonlinear diffusion 
model (|2.2|l and (|2.3l l is independent of the diffusion function b, and the closure of 
S{u-) coincides with the shock set characterized by Oleinik entropy inequalities (or, 
equivalently. Lax shock inequalities). 



Remark 2.3. The conclusions of Theorem \2.2\ do not hold for more general flux- 
functions. This is due to the fact that a strict inequality is required in l|2.91 for the 
existence of the traveling waves. The set based on traveling waves may he strictly 
smaller than the one based on Oleinik entropy inequalities. 



3 Kinetic functions associated with cubic flux-func- 
tions 

Investigating traveling wave solutions of diffusive-dispersive regularizations 
of l|2.Hl is considerably more involved than what was done in Section 2. Besides 
proving the existence of associated (classical and nonclassical) traveling waves 
our main objective will be to derive the corresponding kinetic functions for 
nonclassical shocks. 

To explain the main difficulty and ideas it will be useful to treat first, in the 
present section, the specific diffusive-dispersive model with cubic flux 

dtu + dxV? = e Uxx + 6 Uxxx, (3.1) 

which, formally as e, 6^0, converges to the conservation law with cubic flux 

dtu + dxU^ = 0. (3.2) 

We are interested in the singular limit e ^ in ||3.1| | when the ratio 

a = 4= (3.3) 

yd 

is kept constant. We assume also that the dispersion coefficient 5 is positive. 
Later, in Theorem 14.51 below, we will see that all traveling waves are classical 
when 6 < which motivates us to restrict attention to 6 > 0. 

We search for traveling wave solutions of l|3.1l l depending on the rescaled 

variable 

X- At X- At ,^., 

y:=a^ = -^. (3.4) 

Proceeding along the same lines as those in Section 2 we find that a traveling 
wave y i—> u{y) should satisfy 

- A My + (U\ = aUyy + Uyyy, (3.5) 

together with the boundary conditions 



lim ii{y) = u±, 

lim Uy{y) = lim Uyy{y) = 0, 

y—>±co ■' y—>±co ■'■' 



(3.6) 



where U- + u+ and A are constants. Integrating ( 13.51 1 once we obtain 

a Uyiy) + Uyyiy) = -A {u{y) - uJ) + u{yf -ut, ye R, (3.7) 

which also implies 

H — tl 

A = — =ul + u-u+ + ul. (3.8) 

U+-U- 



To describe the family of traveling waves it is convenient to fix the left-hand 
state (with for definiteness U- > 0) and to use the speed A as a parameter. Given 
U-, there is a range of speeds, 

Ae{3ul/A,3ul), 

for which the line passing through the point with coordinates (m-, ut) and with 
slope A intersects the graph of the flux f{u) :- «■' at three distinct points. For the 
discussion in this section we restrict attention to this situation, which is most 
interesting. There exist three equilibria at which the right-hand side of (|3.7b 
vanishes. The notation 

M2 < Ml < Mo := M- 

will be used, where U2 and Mi are the two distinct roots of the polynomial 

u^ + Mo M + Mq = A. (3.9) 

Observe in passing that U2 + u\ + uq - 0. 

Consider a trajectory y i-^ u{y) leaving from m_ at -co. We want to determine 
which point, among Mi or U2, the trajectory will reach at +oo. Clearly, the 
trajectory is associated with a so-called classical shock if it reaches Mi and with 
a so-called nonclassical shock if it reaches M2. Accordingly, we will refer to it as 
a classical trajectory or as a nonclassical trajectory, respectively. 

We reformulate ll3.7b as a differential system of two equations, 

AQ = K(M,.), (3.10) 

where 

K{u,v) = i , / ^^ / ,^|, g{u,A) = u^-Au. (3.11) 

^ ' \-av + g(u,A) - g(u-,A)j &^ ' ' ^ ' 

The function K vanishes precisely at the three equilibria (mq, 0), (mi, 0), and (m2, 0) 
of | |3.10|| . The eigenvalues of the Jacobian matrix of K{u, v) at any point (u, 0) are 

-a/2 ± ^ja^l^ + g[^{u,A). So we set 

^{u) = \(-a- V«^ + 4(3m2-A)), 

^ (3.12) 



1 

Jl{u) =j(-a+ V«^ + 4(3m2-A)) . 



At this juncture, we recall the following standard definition and result. (See the 
bibliographical notes for references.) 

Definition 3.1 (Nature of equilibrium points). Consider a differential system of the 
form (|3.101 where K is a smooth mapping. Let (m„ c.) e R^ he an equilibrium point, 
that is, a root of K(u,,v,) = 0. Denote by fi = n{u,,v,) and Ji = Ji{u,,v,) the tzvo 
(real or complex) eigenvalues of the Jacobian matrix ofK at (m„ v,), and suppose that a 
basis of corresponding eigenvectors r{u,, v,) and r(u,, v,) exists. Then, the equilibrium 
(m„c,) is called 



1. a stable point ifRe(ji) and Re(Jl) are both negative, 

2. a saddle point ff Re(/,i) and ReQl) have opposite sign, 

3. or an unstable point ifRe{p.) and ReQI) are both positive. 

Moreover, a stable or unstable point is called a node if the eigenvalues are real and a 
spiral if they are complex conjugate. 

Theorem 3.2 (Local behavior of trajectories). Consider the differential system (|3.101 
under the same assumptions as in Definition 2.1. If(u,,v,) is a saddle point, there 
are two trajectories defined on some interval {-co, y,) and two trajectories defined on 
some interval (y,, +oo) and converging to (m„ v,) at -co and +co, respectively. The 
trajectories are tangent to the eigenvectors r{u,, v,) and r{u,, v,), respectively. 

Retttrnrng to (|3.11t and ||3.12|| we conclude that, since ^'„(m. A) = 3u^ - A is 
positive at both u = U2 and u - Uq, we have 

p{Uo) < < JI(Uo), jU(M2) < < Ji{U2). 

Thus both points M2 and mq are saddle points. On the other hand, since we have 
g,',(Mi, A) < 0, the point mi is stable: it is a node if a^ + 4 (3 Mj - A) > or a spiral 
if a^ + 4 (3 Mj - A) < 0. In summary, for the system (|3.101 - l|3.111 

U2 and Un are saddle points and 

^ (3.13) 

Ml is a stable point (either a node or a spiral). 

In the present section we check solely that, in some range of the parameters 
Mo, A, and a, there exists a nonclassical trajectory connecting the two saddle 
points Mo and M2. Saddle-saddle connections are not "generic" and, as we will 
show, arise only when a special relation (the kinetic relation) holds between mo, 
A, and a or, equivalently, between uq, M2, and a; see l|3.15|l below. 

For the cubic model O.ll l an explicit formula is now derived for the nonclas- 
sical trajectory. Motivated by the fact that the function g in B.llll is a cubic, we 
a priori assume that v - Uy is a parabola in the variable u. Since v must vanish 
at the two equilibria we write 

v{y) = a{u{y)-U2)(u{y)-UQ), y e R, (3.14) 

where a is a constant to be determined. Substituting l|3.14b into (|3.101 l- (|3.111 l, we 
obtain an expression of Vy: 

Vy = -av + u^ - Uq - A {u - Mo) 

= -av + {u- M2) (m - mq) (m + Mo + "2) 



= V (-a + - (m + Mo + M2)). 



But, differentiating (|3.14l l directly we have also 

Vy - aUy{2u-uo- M2) 
= av{2u- uq- ux). 

The two expressions of Cy above coincide if we choose 

1 1 

- =2a, -a + - (uq + M2) = -a{uo + U2). 
a a 

So, a = 11 V2 (since clearly we need c < 0) and the three parameters uq, U2, and 
a satisfy the explicit relation 

V2 
M2 = -Mo + -IT- «• (3.15) 

Since mj = -uq - U2 we see that the trajectory (|3.14l l is the saddle-saddle con- 
nection we are looking for, only if 112 < Mi as expected, that is, only if 

Mo > —7— a. (3.16) 

Now, by integrating (|3.14l l, it is not difficult to arrive at the following explicit 
formula for the nonclassical trajectory: 



U0 + U2 U0-U2 ^ ,/ Mo - M2 \ 

. ^ - („_ _ ^L) tanh((M_ - ^ A). 
3V2 ^ 3V2^ ^^ 3V2^ V2^ 



(3.17) 



We conclude that, given any left-hand state Mo > 2 y2a/3, there exists a saddle- 
saddle connection connecting mq to — Mq + y2a/3 which is given by l|3.171 . 
Later, in Section 4 and followings, we will prove that the trajectory just found is 
actually the oniy saddle-saddle trajectory leaving from mo > 2 y2a/3 and that 
no such trajectory exists when mo is below that threshold. 

Now, denote by S(,(m-) the set of all right-hand states u+ attainable through 
a diffusive-dispersive traveling wave of O.ll l with 6 > and e/ v6 = a fixed. 
In the case of the equation l|3.1l l the results to be established in the following 
sections can be summarized as follows. 

Theorem 3.3 (Kinetic function and shock set for the cubic flux). The kinetic 
function associated with the dijfusive-dispersive model B.ll is 



<pl(u.) 



-u- - a/2, u- < -a, 

-U-/2, |m_| < a, (3.18) 

-M_ + a/2, u- > a, 



with a •.= 2 a V2/3, while the corresponding shock set is 
S„(m_) = 



'(m_, a/2] U { - M_ - a/2}, m_ < -a, 
[-u-/2,u-), -a<u-<a, (3.19) 

{ - M_ + a/2} U [-a/2, m_), m_ > a. 



In agreement with the general theory of the kinetic function, | |3.18l l is mono- 
tone decreasing and lies between the limiting functions (p^{u) := -ii/2 and 
(pgiu) := —u. Depending on m_ the shock set can be either an interval or the 
union of a point and an interval. 

Consider next the entropy dissipation associated with the nonclassical shock: 

E(M_; a, U) := - ((p^iu.f + cp^iu.) u. + uVj (!i(<p!,("-)) - U{u.)) 
+ f((p^(w_))-F(w_), 

where {U,F) is any convex entropy pair of the equation | |3.2| |. By multiplying 
ll3.5l l by U'{u{y)) and integrating over y e R we find the equivalent expression 



E{iu;a,U)= i U'{u{y))(aUyy{y) + Uyyy{y))dy 
Jr 

= f(-aU"{u)ul + U"'{u)ul/2)dy. 
Jr 



(3.21) 



So, the sign of the entropy dissipation can also be determined from the explicit 
form | |3.17|| of the traveling wave. 

Theorem 3.4 (Entropy inequalities). 1. For the quadratic entropy 

U{z) = 2^12, z e ]R, 
the entropy dissipation E{u-; a, U) is non-positive for all real U- and all a > 0. 

2. For all convex entropy U the entropy dissipation £(«_; a, U) is non-positive for 
alla>0 and all |m_| < 2 V2a/3. 

3. Consider |m_| > 2 V2a/3 and any (convex) entropy U whose third derivative is 
sufficiently small, specifically 

[\u-\ - a/{3 V2))V"'(z)l < 2a V2 !J"(z), z e M. (3.22) 

Then, the entropy dissipation E(m_; a, U) is also non-positive. 

4. Finally given any \u-\ > 2 V2a/3 there exists infinitely many strictly convex 
entropies for which £(«_; a, U) is positive. 



Proof. When U is quadratic (with U" > and W" = 0) we already observed 
that Item 1 follows immediately from (|3.21l l. The statement Item 2 is also 
obvious since the function cp'' reduces to a classical value in the range under 
consideration. Under the condition | |3.22t the integrand of l|3.21t is non-positive, 
as follows from the inequality (see l|3.14t ) 

\uy\ <^{uq- mf = ^(u-- a/(3 V2)) . 
4V2 V2^ ' 

This implies the statement Item 3. Finally, to derive Item 4 we use the (Lipschitz 
continuous) Kruzkov entropy pairs 

Uk{z) ■- \z - k\, Ffc(z) := sgn(z - k){r' - k^), z e IR, (3.23) 

with the choice k = -m_/2. We obtain 

E(u-;a, Uk) = ^ |m-| (|m-| - 2a V2/3)^ > 0. 

By continuity, E(m_; a, Lf/c) is also strictly positive for all fc in a small neigh- 
borhood of -M_/2. The desired conclusion follows by observing that any 
smooth convex function can be represented by a weighted sum of Kruzkov 
entropies. D 

Remark 3.5. We collect here the explicit expressions of some functions associated with 
the model (|3.11 l. From now on we restrict attention to the entropy pair 

U{u) = u^l2, F{u) = 3u^/A. 

First of all, recall that for the equation l|3.2b the following two functions 

cp\u) = -'^, cpl{u) = -u, weR. (3.24) 

determine the admissible range of the kinetic functions. 
PVe define the critical diffusion-dispersion ratio 

3 
A{uq, U2) = — = (mo + "2) (3.25) 

V2 

for uq > and 1/2 e (-mq, -mo/2) and for uq < and M2 e (-mo/2, -mq)- In view 
of Theorem 13.31 (see also (|3.151 ), a nonclassical trajectory connecting uq to 1/2 exists 
if and only if the parameter a = ej ^fb equals A{uq,U2). The function A increases 
monotonically in M2 from the value to the threshold diffusion-dispersion ratio 
(uo > 0) 

A\uo) = ^. (3.26) 

2 V2 

For each fixed state mq > there exists a nonclassical trajectory leaving from uq if 

and only if a is less than A^(mo)- On the other hand, for each fixed a there exists a 

nonclassical trajectory leaving from wq if and only if the left-hand state mq is greater 

than A^ {a). The function J^ is a linear function (for uq > 0) with range extending 

therefore from A^ = to A^ = +00. 



10 



Remark 3.6. It is straightforward to check that if (|3.1t is replaced with the more 
general equation 

dtu + dJK M^) = e W.VX + 6C M„x, (3.27) 

where C and K are positive constants, then l |3.26l l becomes 

A\uo) = ^ylKC. (3.28) 

2 V2 

Remark 3.7. Clearly, there is a one-parameter family of traveling waves connecting 
the same end states: If u = u{y) is a solution of (|3.5b and Il3.6t , then the translated 
function u = u{y + b) (b e R) satisfies the same conditions. However, one could show 
that the trajectory in the phase plane connecting two given end states is unique. 

4 Kinetic functions associated with general flux-func- 
tions 

Consider now the general diffusive-dispersive conservation law 

dtu + drf{u) = £ {b{u) Mx)^. + 6 (ci(m) (c2(m) u^)x)^, u = u'-\x, t), (4.1) 

where the diffusion coefficient b{u) > and dispersion coefficients c\{u), C2{u) > 
are given smooth functions. We assume that / : ]R ^ ]R is a concave-convex 
function satisfying, by definition. 



uf"{u)>0 for all M 7^0, 
f"'{0)*0, lim f'{u) = +oo. 



(4.2) 



We are interested in the singular limit e ^ when 6 > and the ratio a = e/ yd 
is kept constant. The limiting equation associated with l|4.1t , formally, is the 
scalar conservation law 

dfU + dxfiu) = 0, u = u{x, t) e ]R. 

It can be checked that the entropy inequality 

dtU{u) + d,F{u) < 

holds, provided the entropy pair (U, F) is chosen such that 

U"{u) := ^, F'(m) := U'{u)f'{u), m e M, (4.3) 

which we assume in the rest of this paper. Since ci,C2 > the function LI is 
strictly convex. 
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Given two states u+ and the corresponding propagation speed 

A = a{u-, u+) 



'fJ^^^tns^ u,^u., 



U+-U- 



yf'iu-), u+ = U-, 

we search for traveling wave solutions u - m(i/) of l|4.1l l depending on the 
rescaled variable y := (x - A t) a/e. Following the same lines as those in Sec- 
tions 1 and 2 we find that the trajectory satisfies 

Ci(u){c2{u)Uij)y + ab{u)Uij = -/\{u-u-) + f{u) - f{u-), u = u{y), (4.4) 

and the boundary conditions 

lim u{y) = u±, lim Uy{y) = 0. 

Setting now 

V = C2{U) Uy, 

we rewrite l|4.4) l in the general form (|3.10t for the unknowns u = u{y) and 
^ = ^(y) (y e R), i.e., 

!(:)=«„,., (4.5) 



with 

X(M, ^) - _ b(") ,, , g(»,A)-g(»-,A) 

V " Cl(u)C2(») '^ "^ Ci(») 

while the boundary conditions take the form 



g{u,A):=f{u)-\u, (4.6) 



lim u{y) = u±, lim v{y) = 0. (4.7) 

y— *±oa y— >±oo 

The function K in ( |4.6l l vanishes at the equilibrium points {u, v) e R^ satis- 
fying 

g{u,A) = g{u.,A), v = 0. (4.8) 

In view of the assumption (|4.2t , given a left-hand state «_ and a speed A there 
exist at most three equilibria u satisfying l|4.81 (including m_ itself). Considering 
a trajectory leaving from m_ at -oo, we will determine whether this trajectory 
diverges to infinity or else which equilibria (if there is more than one equilibria) 
it actually connects to at +oo. Before stating our main result (cf. Theorem I4.3[ 
below) let us derive some fundamental inequalities satisfied by states U- and 
u+ connected by a traveling wave. 
Consider the entropy dissipation 

E{u-, M+) ■- -a{u-, u+) (U{u+) - U(m_)) + F(m+) - F(m_) (4.9) 
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or, equivalently, using (|4.31 and (I4.7I > 

X+oo 
!J'("(y)) (-A My(y) + f{u(y))y) dy 
CO 

X+co 
U"{u{y)) (-A (M(y) - w_) + /(w) - /(«_)) My(y) dy (4.10) 
CXI 






,fl(M_,M+)) -^(M_,fl(M_,U+))) ^— dz. 



Ci(z) 

In view of 

X+CX3 
!J'(m) (a (b{u) My)^ + (ci(m) (c2(m) My)y)^ ) rfy 

aU"{u)h{u)u]dy, 

we have immediately the following. 

Lemma 4.1 (Entropy inequality). If there exists a traveling wave o/l l4.4|l connecting 
U- to u+, then the corresponding entropy dissipation is non-positive, 

E{u-,u+) < E(u-,U-) = 0. 

From^ the graph of the function / we define the functions (p^ and A'' by 

^ ' u - (p\u) 

We have u <p^(u) < and by continuity (p^{0) = and, thanks to (|4.2| |, the map 
(p^ : ]R ^ ]R is decreasing and onto. It is invertible and its inverse function is 
denoted by tp"*". Observe in passing that, m_ being kept fixed, A''(m-) is a lower 
bound for all shock speeds A satisfying the Rankine-Hugoniot relation 

-A(m+-m-) + /(m+)-/(m-) = 

for some u+. 

The properties of the entropy dissipation l|4.9|l are determined from the 



zero-entropy dissipation function (pj! was introduced. 



Lemma 4.2 (Entropy dissipation function). There exists a decreasing function 
cpp : M ^ M such that for all U- > (for instance) 

E{u-, u+) = and u+ + m_ if and only if u+ = (pg(M_), 

E(m_,m+) < if and only if (fgiu-) < u+ < m_, 

and 

cp-\u.) < cpliu.) < (p\u.). 
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In passing, define also the function cp' = (Pg(M-) and the speed Aq = Ao(m-) 
by 

/(«-)- M(u-)) f{u.)-f{4(u.)) 

'^o^"-) = ^r^ — = i ' "- ^ 0- (4-11) 

Combining Lemmas 14.11 and 14.21 together we conclude that, if there exists a 
traveling wave connecting m_ to m+, necessarily 

u+ belongs to the interval [^q(m_), m_]. (4-12) 

In particular, the states u+ > U- and u+ < (p~''(M_) cannot be reached by a 
traveling wave and, therefore, it is not restrictive to focus on the case that three 
equilibria exist. 

Next, for each m_ > we define the shock set generated by the diffusive- 
dispersive model l|4.Hl by 

§a{u-) := |m+ / there exists a traveling wave of l|4.41 connecting m_ to m+I. 

Theorem 4.3 (Kinetic function and shock set for general flux). Given a concave- 
convex flux-function f (see H4.2I I), consider the diffusive-dispersive model (|4.1| | in which 
the ratio a = e/ yd > is fixed. Then, there exists a locally Lipschitz continuous and 
decreasing kinetic function (p\ : IR ^ IR satisfying 



(4.13) 



(p\u) < <(w) < <p»(m), m < 0, 
(Po(") < <?«(") ^ (p\u), M > 0, 
and such that 

Here, the function (p„ is defined from the kinetic function cp), by 
f{u) - ficpiiu)) f{u)-f(<pi{u)) 
u - (pl{u) u - (pl{u) 

with the constraint 

cpliu) < cpliu) < (p\u), M < 0, 

(p\u) < <pi{u) < <pI{u), m > 0. 
Moreover, there exists a function 

A^:U^ [0, +co), 

called the threshold diffusion-dispersion ratio, which is smooth away from u = 0, 
Lipschitz continuous at u = 0, increasing in m > 0, and decreasing inu <0 with 

A^(m)~C|w| flSM^O, (4.16) 
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(4.15) 



(where C > depends upon f, b, cj, and C2 only) and such that 

(p^aiu) = (p\u) when a > A\u). (4.17) 

Additionally we have 

(p\{u) -^ (pl(u) as a-^0 for each u e R. (4.18) 



The proof of Theorem 14.31 will be the subject of Sections 4 and 5 below. 
The kinetic function (p|^ : ]R ^ ]R completely characterizes the dynamics of the 
nonclassical shock waves associated with ||4.1| |. In view of Theorem l4.3l one can 
solve the Riemann problem. The kinetic function (p„ is decreasing and its range 
is limited by the functions q)^ and (p„. Therefore we can solve the Riemann 
problem, uniquely in the class of nonclassical entropy solutions selected by the 
kinetic function (p[', . 

The statements (|4.171 and (|4.181 provide us with important qualitative prop- 
erties of the nonclassical shocks: 

1. The shocks leaving from m_ are always classical if the ratio a is chosen to 
be sufficiently large or if m_ is sufficiently small. 

2. The shocks leaving from U- are always nonclassical if the ratio a is chosen 
to be sufficiently small. 

Furthermore, under a mild assumption on the growth of / at infinity, one could 
also establish that the shock leaving from m_ are always nonclassical if the state 
U- is sufficiently large. (See the bibliographical notes.) 

In this rest of this section we introduce some important notation and in- 
vestigate the limiting case when the diffusion is identically zero (a = 0). We 
always suppose that m- > (for definiteness) and we set 

Mo = U-. 

The shock speed A is regarded as a parameter allowing us to describe the set of 
attainable right-hand states. Precisely, given a speed in the interval 

Ae(A\uo),f'{uo)), 

there exist exactly three distinct solutions denoted by uq, mi, and U2 of the 
equation ( |4.8l l with 

M2 < (p\uq) <Ui < Uq. (4.19) 

Recall that no trajectory exists when A is chosen outside the interval limited by 

A^(Mo)and/'(Mo)- 

From Lemmas 14.11 and 14.21 (see l|4.12|l ') it follows that a trajectory either is 
classical if Uo is connected to 

Ml e [<p\uo), Mo] with A e [a^(mo),/'(mo)] (4.20) 
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or else is nondassical if mq is connected to 

"2 e [<p[5(mo), (p^("o)) with A e (a^(mo), Ao(uo)]- (4.21) 

For the sake of completeness we cover here both cases of positive and 
negative dispersions. For the statements in Lemma l44l and Theorem |4.5| below 
only we will set a :- e/ ^f\6\ and rj - sgn(6) = ±1. If {u,v) is an equilibrium 
point, the eigenvalues of the Jacobian matrix of the function K{u, v) in (|4.6t are 
found to be 




b{u) - b{uf ^ f'{u) - A 



2\ Ci(m)C2(m) \ Ci(m)2c2(m)2 Cl{u) C2{u) 

So, we set 



r^ b(M) / _ . _ 1 , l!l ci{u)c2{ u) 
2 ci(m)c2(m)\ 'V «' Huf 



f^(«; '^' «) = V ^7TX^7T:^ I - 1 -'] a/1 + :i ^T^f^ (/'(")- ^^) 



(4.22) 



J|a Mm) / ^ 1+1^ ci(»)c2( m) 

2 Ci(m)C2(m) \ '' V a2 ^(„)2 



f^(«; A,a) = ^ ::-£^^ I - 1 + '] a/1 O ^^^^S^ (/'(«) - '^) 



Lemma 4.4 (Nature of equilibrium points). Fix some values U- and A and denote 
by (u„ 0) any one of the three equilibrium points satisfying (|4.8|l . 

1. Ifrj = +1 and f'{u,) - A < 0, then (m„ 0) z's a stable point. 

2.1frj {f'{u,) - A) > 0, then (m„ 0) z's a saddle point. 

3. Ifri = -1 and f'{u,) - A > 0, then {u„ 0) z's an unstable point. 

Furthermore, in the two cases that ri{f'{u,) - A) < we have the additional result: 
When a^ b{u,)^ + 4 t] ci{u,) C2{u,) {f'{u,) - A) > the equilibrium is a node, and z's a 
spiral otherwise. 

For negative dispersion coefficient b, that is, when i] - -1, we see that both 
Ui and M2 are unstable points which no trajectory can attain at +oo, while ui is 
a stable point. So, in this case, we obtain immediately: 

Theorem 4.5 (Traveling waves for negative dispersion). Consider the diffusive- 
dispersive model \A.\\ where the flux satisfies \A.2\ . If e > and 6 < 0, then only 
classical trajectories exist. 

Some additional analysis (along similar lines) would be necessary to estab- 
lish the existence of these classical trajectories and conclude that 



§„(m_) = §(m_) := 



'[(pH"-),"-], M->0 

when 6 < 0, 

[lZ_,(p''(M_)], M_ <0 
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which is the shock set already found in Section 2 when 6 = 0. 

We return to the case of a positive dispersion which is of main interest here. 
(From now on 7] = +1.) Since g[i{^, A) is positive at both u = U2 and u = Uq, we 
have 

|U(Mo) < < JI(Uo), /j(M2) < < f/(M2), 

and both points W2 and Uq are saddle. On the other hand, since g'„(wi. A) < 0, 
the equilibrium Hi is a stable point which may be a node or a spiral. These 
properties are the same as the ones already established for the equation with 
cubic flux. The following result is easily checked from the expressions l|4.22|l . 

Lemma 4.6 (Monotonicity properties of eigenvalues). In the range of parameters 
where n{u, A, a) and Jl{u; A, a) remain real-valued, we have 

d^ d^ gjl 

-=rr(u; A, a) > 0, r—{u; A, a) < r—(u; A, a) < 0, 
c/A da oA 

and, under the assumption f'{u) - A > 0, 

dTi 

-T^{u;A,a) < 0. 
da 

To the state uq and the speed A e (a''(mo)/ Ao(mo)) we associate the following 
function of the variable u, which will play an important role throughout, 

rV \ C2(2) 

G(m; Uq, A) := [g{z. A) - ^(mo, A)j — — dz. 

Observe, using I I4.10K that the functions G and E are closely related: 

G{u; uo. A) = -E{uq, u) when A = a(uo, u). (4.23) 

Note also that the derivative di,G{u; uq. A) vanishes exactly at the equilibria uq, 
u\, and M2 satisfying ||4.8|I . Using the function G we rewrite now the main 
equations I l4.51 - (|4.6l l in the form 

C2(w) Uy = V, (4.24) 

C2iu)Vy = -a-y{v + G'^iu;uo,A), (4.25) 

Ci(U) 

which we will often use in the rest of the discussion. 

We collect now some fundamental properties of the function G. 

Theorem 4.7 (Monotonicity properties of the function G). Fix some uq > and 
A e (a''(mo)//'(wo)) ^f^d consider the associated states Ui and U2. Then, the function 
u 1-^ G{u) := G{u; ug. A) satisfies the monotonicity properties 

G'{u) < 0, u < U2 or u e (wi, uq), 
G'{u) > 0, u e (m2. Ml) or u > uq. 

17 



Moreover, if A e ( A''(mo)/ Ao(mo)) ^^ ^^^^ 

G(mo) = < G(U2) < G(ui), (4.26) 

w/zf/e, r/A = Ao(uo)/ 

G(mo) = G(m2) = < G(mi) (4.27) 

and finally, if A e (Ao(mo)//'(wo))/ 

G(m2) < = G(uo) < G(ui). (4.28) 

Proof. The sign of G' is the same as the sign of the function 

g{u. A) - g{uo. A) = {u- Mo) ( - a). 

^ U — Mo 

So, the sign of G' is easy determined geometrically from the graph of the 
function /. To derive | |4.26| |-| |4.28| | note that G(mo) - and (by the monotonicity 
properties above) G(mi) > G(uo)- To complete the argument we only need the 
sign of G(u2)- But by l|4.231 we have G(m2) - -E{uo, 1/2) whose sign is given by 
Lemma 10 D 

We conclude this section with the special case that the diffusion is zero. 
Note that the shock set below is not the obvious limit from | |4.14| |. 

Theorem 4.8 (Dispersive traveling waves). Consider the traveling wave equation 
l|4.4|l in the limiting case a = (not included in Theorem l4.3l > under the assumption 
that the flux f satisfies (|4.2)l . Then, the corresponding shock set reduces to 

§o(m-) = j(p[;(M-),u-}, u-eR. 

Proof. Suppose that there exists a trajectory connecting a state m_ > to a state 
M+ i= U- for the speed A = a{u-, u+) and satisfying (see (|4.241 - (|4.25t ) 

Cliu) Uy = v, 

■' (4.29) 

ci{u)Vy = g{u,A)-g{u-,A). 

Multiplying the second equation in l|4.29t by vlci{u) = C2{u) Uy/c\{u), we find 

and, after integration over some interval (-00, y], 

^vHy) = G{u{yy,u.,A), yeU. (4.30) 

Letting y -^ +(x> in | |4.30|| and using that v(y) -^ we obtain 

G(m+;m-,A) = 
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which, by l|4.231 , is equivalent to 

E(m+, M_) = 0. 



Using Lemma [4. 21 we conclude that the right-hand state u+ is uniquely deter- 
mined, by the zero-entropy dissipation function: 

M+ = <p|5(u_), A = Ao(m-). (4.31) 

Then, by assuming (|4.311 and U- > 0, Theorem l4.7l implies that the function 
u H-* G(m; m_. A) remains strictly positive for all u (strictly) between u+ and U-. 
Since v <0 we get from | |4.30| | 

v{y) = - V2G(M(y);M_,A). (4.32) 

In other words, we obtain the trajectory in the (m, v) plane: 

V = v{u) = - y2G(M7Ml7A), u e [u+, m_], 
supplemented with the boundary conditions 

v{u-) = v{u+) = 0. 

Clearly, the function v is well-defined and satisfies v{u) < for all u e (m+, m_). 
Finally, based on the change of variable y e [-oo, +oo] h^ m = ^(y) e [u^, «_] 
given by 

dy = _, , dw, 
i?(m) 

we immediately recover from the curve v = v{u) the (unique) trajectory 

yi-^("(y)/%))- 

This completes the proof of Theorem |4.8l n 

5 Traveling waves corresponding to a given speed 

We prove in this section that, given Uq, U2, and A - a{uQ, 1/2) in the range (see 

623) 

U2 e [cpliuo), (p\uo)), A e (a\uo), Ao(mo)], (5.1) 

a nonclassical connection always exists if the ratio a is chosen appropriately. 
As we will show in the next section this result is the key step in the proof of 
Theorem 14.31 The main existence result proven in the present section is stated 
as follows. 
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Theorem 5.1 (Nonclassical trajectories for a fixed speed). Consider two states 
uq > and ux <0 associated with a speed 



A = a(uo,U2) e (aVo)/ Ao(mo)]- 



Then, there exists a unique value a >0 such that Uq is connected to 112 by a diffusive- 
dispersive traveling wave solution. 



By Lemma 1441 uq is a saddle point and we have Ji{uo) > and from Theo- 
rem l3.2l it follows that there are two trajectories leaving from Uq aty - -00, both 
of them satisfying 

v{y) _ 

lim — ^ p{uo;A,a)c2iuo). (5.2) 

y^-oo u{y) - Uq 

One trajectory approaches (mq/O) in the quadrant Qi = |w > Uq, v > Ol, the 

other in the quadrant Q2 = \u < uq, v < Ol. On the other hand, M2 is also a 

saddle point and there exist two trajectories reaching M2 at y = +00, both of 

them satisfying 

v{y) 

lim —r- = p{u2;A,a)c2iu2). (5.3) 

y^+oo u{y) -U2 — 

One trajectory approaches {u2, 0) in the quadrant Q3 = Im > M2/ ^ < Ol, the other 
in the quadrant Qt = \u < U2, v > 0|. 

Lemma 5.2. A traveling wave solution connecting uq to M2 must leave the equilibrium 
(mq, 0)aty = -co in the quadrant Q2, and reach {u2, 0) in the quadrant Q^aty = +00. 

Proof. Consider the trajectory leaving from the quadrant Qi, that is, satisfying 
u > Uq and c > in a neighborhood of the point (mq/O). By contradiction, 
suppose it would reach the state M2 at +00. Since M2 < Uq by continuity there 
would exist yo such that 

u{yo) = Mo- 
Multiplying l|4.25t by My = v/c2(u) we find 

(v^/2) + a .y , . v^ = G;(m; mq, A) Uy. 
^ 'y ci{u)C2(u) ^ 

Integrating over (-00, yo] we arrive at 

^Hyo) r" 7 b(u) 

Therefore v{yo) = and, since M(yo) - mq, a standard uniqueness theorem for 
the Cauchy problem associated with (|4.241 - (|4.251 l implies that u = uq and v = 
on R. This contradicts the assumption that the trajectory would connect to M2 
at +00. 
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The argument around the equihbrium {112, 0) is somewhat different. Suppose 
that the trajectory satisfies M < M2andi' > in a neighborhood of the point (!i2,0). 
There would exist some value 1/1 achieving a 7oca7 minimwri, that is, such that 

u(i/i) < U2, Uyiyi) = 0, Uyyiyi) > 0. 

From (|4.24) l we would obtain u(i/i) = and, by differentiation of l|4.24|l , 

5^y(l/l) = Wyy(yi)C2(M(l/l)) > 0. 

Combining the last two relations with (14.251 1 we would obtain 

G;,(m(i/i);mo,A) > 



which is in contradiction with Theorem|4]7lsinceM(i/i) < U2cLndG[,{u{yi);uo,A) < 0. 

D 

Next, we determine some intervals in which the traveling waves are always 
monotone. 

Lemma 5.3. Consider a trajectory u = u{y) leaving from uq at -oo and denote by ^ 
the largest value such that ui < u{y) < uofor all y e (-00, E,) and u{^) = ui. Then, 
we have 

Uy <0 on the interval (-00, Q. 

Similarly, ifu = u{y) is a trajectory connecting to U2 at +00, denote by £, the smallest 
value such that M2 < u{y) < u\for all y e {£,, +00) and u{E.) = u\. Then, we have 

Uy < on the interval {£,, +00). 

In other words, a trajectory cannot change its monotonicity before reaching 
the value mi. 

Proof. We only check the first statement, the proof of the second one being 
similar. By contradiction, there would exist 1/1 e {-co,E,) such that 

Uy(yi) = 0, Uyyiyi) > 0, Ml < m(i/i) < Mo- 
Then, using the equation l|4.25t would yield G[,{u{yi);uo,A) > 0, which is in 
contradiction with the monotonicity properties in Theorem 14. 71 n 

Proof of Theorem \5J\ For each a > we consider the orbit leaving from Mq and 
satisfying m < Mq and i; < in a neighborhood of (mq, 0). This trajectory reaches 
the line \u = uA for the "first time" at some point denoted by (mi, V-{a)). In 
view of Lemma [531 this part of trajectory is the graph of a function 

[mi,mo] 3 u h^ V-{u;A,a) 

with of course V-{u\;A,a) - V-{a). Moreover, by standard theorems on 
differential equations, V- is a smooth function with respect to its argument 

(m; a, a) e [mi, mq] x (a^(mo), Ao(mo)] x [0, +00). 
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Similarly, for each a > we consider the orbit arriving at M2 and satisfying 
u > U2 and Z7 < in a neighborhood of {U2,0). This trajectory reaches the line 
ju = uil for the "first time" as y decreases from +00 at some point (mi, V+{a)). 
By Lemma [5. 31 this trajectory is the graph of a function 

[m2. Ml] 3 u h^ v+{u;A,a). 

The mapping v+ depends smoothly upon {u,A,a) e [w2/Wi] x (A^(mo)/Ao(mo) X 
[0, +00). 

For each of these curves u h-> V-{u) and u h-> v+{u) we derive easily from 
||4.24|| -| |4.25|| a differential equation in the (m, v) plane: 

v{u) —{u) + a — — - v{u) = g;(m, uq, A). (5.5) 

au ci(m) 

Clearly, the function 

a e [0, +00) 1-^ W(a) : = v+{ui; A, a) - v-{ui; A, a) 
= V4a) - y_(a) 

measures the distance (in the phase plane) between the two trajectories at 
u = Ml. Therefore, the condition W{a) = characterizes the traveling wave 
solution of interest connecting mq to M2. The existence of a root for the function 
W is obtained as follows. 

Case 1: Take first a = 0. 

Integrating I l5.5b with v = V- over the interval [mi, Mq] yields 

1 , 

-{V-{Q)Y = G(mi;mo,A) - G{uo;no,A) = G(mi;mo,A), 

while integrating | |5.5t with v = v+ over the interval [M2, Mi] gives 
-(V+(0))2 = G(mi;mo,A) - G(m2;mo, A). 

When A i= Ao(Mo)/SinceG(M2;Mo, A) > (Theorem l4.7t and V^(a) < (Lemma |5.3l l 
we conclude that W(0) > 0. When A = Ao(mo) we have G(h2; mq. A) = and 
W(0) = 0. 

Case 2: Consider next the limit a -^ +C0. 

On one hand, since V- < 0, for a > we get in the same way as in Case 1 

i(y_(a))2<G(Mi;Mo,A). (5.6) 

On the other hand, dividing ||5.5|I by c = J7+ and integrating over the interval 
[U2/Mi] we find 

v f , r ^(") ^ ^ r"' G;(M;Mo,A) ^ 

y+(a) = -a ——du+ — — du. 

Ju2 ci(m) J„^ v+{u) 
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Since v = C2{u) Uy <0 and G[,(m) > in the interval [u2, mi] we obtain 

V+(a) < -K a (Ml - U2), (5.7) 

where k = inf„g[„2 „j] b{u)/ci{u) > 0. Combining | |5.6| | and I |5.7| | and choosing a 
to be sufficiently large, we conclude that 

W(a) = V+(a) - V-(a) < 0. 

Hence, by the intermediate value theorem there exists at least one value a 
such that 

W{a) = 0, 

which establishes the existence of a trajectory connecting uq to M2- Thanks to 
Lemma F5.3l it satisfies My < globally. 

The uniqueness of the solution is established as follows. Suppose that there 
would exist two orbits v = v(u) and v* - v*{u) associated with distinct values a 
and a* > a, respectively. Then, Lemma [4.61 would imply that 

/I(mo; a, a) < /I(mo; A, a), ju{u2; A,a) < jj{u2; A, a). 

So, there would exist M3 e (m2, mq) satisfying 

dv* dv 

v{u3) = v'ius), -j^im) > ^("3). 



Comparing the equations i5.5i satisfied by both v and v*, we get 

^("^)(|("3)-5;(.3)) = (a-a);^.(«3). (5.8) 



Now, since ^(ms) 9^ (the connection with the third critical point (mi,0) is 
impossible) we obtain a contradiction, as the two sides of (|5.8|l have opposite 
signs. This completes the proof of Theorem 15. II n 



Remark 5.4. It is not difficult to see also that, in the proof of Theorem \5m 

a 1-^ V-{a) is non-decreasing (5.9) 

and 

a 1-^ V+{a) is decreasing. (5.10) 

In particular, the function W{a) := V+{a) - V-{a) is decreasing. 

Theorem 5.5 (Threshold function associated with nonclassical shocks). Consider 
the function A = A{uo, M2) which is the unique value a for which there is a nonclassical 
traveling wave connecting uq to M2 (Theorem \5.1h It is defined for uq > and M2 < 
with M2 e (pg(Mo), (p''(Mo)) or, equivalently, Mq e \(pg{u2),(p~''^{u2)j. Then we have the 
following two properties: 
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1. The function A(uo, u^A is increasing in M2 ^nd maps [(Pq(mo)/ <p''("o)) onto some 
interval of the form 0,A''(mo)) where A^uq) e (0, +00]. 

2. The function A is also increasing in mq and maps the interval [<Pq{u2), <p~'^{u2)) 
onto the interval [o, A''(^"^(m2)))- 

Later (in Section 6) the function A will also determine the range in which 
classical shocks exist. From now on, we refer to the function A as the critical 
diffusion-dispersion ratio. On the other hand, the value A^{uq) is called the 
threshold diffusion-dispersion ratio at uq. Nonclassical trajectories leaving 
from Mo exist if and only if a < A^{uo). 

Observe that, in Theorem 15 .51 we have A{uq,U2) -^ when 112 -^ q]g{uo), 
which is exactly the desired property l|4.18b in Theorem l4.3l 

Proof. We will only prove the first statement, the proof of the second one being 
completely similar. Fix mq > and u*2 < 112 < mq so that 

AV,)<A.^'"''-^'°°'<A-.^'"'.'-^'"°'.A,W 

M2 - Mo "2 ~ "0 

Proceeding by contradiction we assume that 

a* := A{uq, u*2) > a := A{uo, U2). 

Then, Lemma|4l6]implies 

Jl{uo; A, a) > Jl{uo; A,a')> Ji{uq; A*, a*). 

Let V = v{u) and v* = v*{u) be the solutions of (|5.51 associated with a and a*, 
respectively, and connecting uq to U2, and uq to u^, respectively. Since u'^ < U2, 
by continuity there must exist some state uj, e {u2, mq) such that 

/ N *, s dv\ ^ dv , ^ 
v{u3) = V (1/3), -j^im) > ^("3)- 

On the other hand, in view of | |5.5| | which is satisfied by both v and v* we obtain 

.(dv\ dv X , fc("3) ,,. ,,, ,C2(M3) 

which leads to a contradiction since the left-hand side is non-positive and the 
right-hand side is positive. This completes the proof of Theorem l5.5l D 

We complete this section with some important asymptotic properties (which 
will establish l l4T6ll -(|4l7 |l in Theoreml43]|. 

Theorem 5.6. The threshold diffusion-dispersion ratio satisfies the following two 
properties: 
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1. A^{uo) < +coforalluQ. 

2. There exists a traveling wave connecting uq to U2 = <P^{uo)for the value a = 
A\uo). 

Proof. Fix Uq > 0. According to Theorem 15.11 given A e (A''(mo)/Ao(mo)] there 
exists a nonclassical trajectory, denoted by m i-^ v{u), connecting uq to some M2 
with 

, /(W2)-/(Mo) ^ h/ X ., s .riix 

A = , M2 < (P^(Mo)/ « = A(M0/"2)- (5-11) 

M2 - Mo 

On the other hand, choosing any state m^ > mq and setting 

,. /K)-/K) . ^, . 
'^ = u* - »' — ' "i = r("o), 
"0 "1 

it is easy to check from | |4.22|| that, for all a* sufficiently large, ^i(m^;A', a*) 
remains real with 

^i(u;;A*,a*)<0. 

Then, consider the trajectory u 1-^ v*{u) arriving at u*^ and satisfying 
lim - «(Mi;A ,a )c2(u.) < 0. 

"^"i M - M* — ^ 'VI' 

i,>i,i 1 

Two different situations should be distinguished. 

Case 1 : The curve v* - v*{u) crosses the curve v = v{u) at some point M3 where 

M*j < M3 < Mo, '^(Ms) = ^'(Ms), 1-("3) > -l-("3)- 

AM aw 

Using the equation | |5.5| I satisfied by both v and u* we get 

(dv* , ^ dv , ^\ , , , b(uj) 
^-(.3)--(M3)) + (a-a)^^ 

<0. 



^("3) (-i-("3) - l-("3)) + («* - a) —A: i'Ims) = G;,(m3; m*. A*) - G;,(m3; mo. A) 



In view of our assumptions, since i^(m3) < we conclude that a < a* in this first 



case. 



Case 2 : v* = v*{u) does not cross the curve v = v{u) on the interval (m*, mo). 

Then, the trajectory v* crosses the M-axis at some point M4 e (u^,mo]- Inte- 
grating the equation ll5.5l l for the function v on the interval [M2, mo] we obtain 

a — -— v{u) du - G(m2; Mo, A) - G(mo; mo. A). 

Juo ci(m) 

On the other hand, integrating (|5.5b for the solution v* over [m* , 1/4] we get 
a' f ' -^ v'{u) du = G(m;; m*. A*) - G(m4; m*. A*). 
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Since, by our assumption in this second case, 

r'" b(u) r"i b(u) 

Juo Ci(m) Ju, Ci(m) 

we deduce from the former two equations that 

, G(m2; mq. A) - Gjup; Up, A) 
'^ - "^ G(w;; M*, A*) - G(m4; Mq' '^*) " ' 

where C is a constant independent of M2- More precisely, M2 describes a small 
neighborhood of q}^{uo), while u* u* M4, and A* remain fixed. 
Finally, we conclude that in both cases 

A(mo, "2) < C a*, 

where a* is sufficiently large (the condition depends on uq only) and C is 
independent of the right-hand state U2 under consideration. Hence, we have 
obtained an upper bound for the function H2 1-^ A{uo,U2)- This completes the 
proof of the first statement in the theorem. 

The second statement is a consequence of the fact that A{uq, U2) remains 
bounded as 1/2 tends to (p^iuo) and of the continuity of the traveling wave v with 
respect to the parameters A and a, i.e., with obvious notation 

v{.;A\uo),A^{uo))= lim v{.;A{uo,U2),A{uo,U2)). 



The function A^ = A''(mo) maps the interval (0, +00) onto some interval 
[A^, A^i] where < A^ < A^ < +00. The values A^ and A^ correspond to 7ower 
and upper bounds for the threshold ratio, respectively. The following theorem 

shows that the range of the function A^ug), in fact, has the form 0, A^i . 

Theorem 5.7. With the notation in Theorem [5^ the asymptotic behavior o/A^(uo) as 
Mo — > is given by 

A\uo)~KUo, k:=^1^|PV37^>0. (5.12) 

Note that of course (|4.2b implies that /'"(O) > 0. In particular. Theorem 15. 71 
shows that A^(0) = A^(0) = 0. Theorem 15.71 is the only instance where the 
assumption /'"(O) # (see jA.l} ) is needed. In fact, if this assumption is 
dropped one still have A^up) ^ as mq — » 0. (See the bibliographical notes.) 

Proof. To estimate A'' near the origin we compare it with the corresponding 
critical fvmction A; determined explicitly from the third-order Taylor expansion 
/* of / = f{u) at u = 0. (See (|5.171 below.) We rely on the results in Section 3, 
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especially the formula B.26I I which provides the threshold ratio explicitly for 
the cubic flux. 

Fix some value uq > and the speed A - A\uo) so that, with the notation 
introduced earlier, U2 = ui = <p\uo). Since /'"(O) 9^ it is not difficult to see that 

U2 = cp\uo) = -(1 + O{uo)) ^ 

(as is the case for the cubic flux f{u) - u^). A straightforward Taylor expansion 
for the function 

G{u) := G{u;uo,A^{uo)) 

yields 

G(m) - G(m2) = G{u) - G{(p\uo)) 



(/'"(O) ^ (3 U2 + U) + 0(|M2P + Iwft). 



24 y ' ' ci(0) 
Since, for all u e [u2, uq] 

4:U2 < U + 3 U2 < Uq + 3 U2 = M2 (1 + 0{Uo)), 

we arrive at 



G(«)-G(«2)-r(0)^(« + 3«2)^"""^^' 



< C Mo I" + 3 M2I (m - W2) ■ (5.13) 



ci(0) ^ '' 24 

Now, given e > 0, we can assume that Mq is sufficiently small so that 

(0 -|(l + ^)<U2<-f(l-.), 

ci(0) ci(m) ci(0) 

(n'O cy(0) (1 - £) < c/m) < Cj{0) (1 + £), u e [m2. Mo], 7 = 1, 2. 

Introduce next the flux-function 

M^ 

Uu) = k—, fc = (l + e)/"'(0), mgR. (5.15) 

6 

Define the following (constant) functions 

V{u) = h{Q), c\{u) = cm, c*(m) = C2(0). 

To these fimctions we can associate a function G, by the general definition in 
Section 4. We are interested in traveling waves associated with the functions 
/,, h* , c*, and c*^, and connecting the left-hand state Mp given by 

iL = -2 M2 
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to the right-hand state H2 (which will also correspond to the traveling wave 
associated with /). 

The corresponding function 

G.(m) := G,{u;Uq,A^{Uq)) 
satisfies 

G.{u) - G.(U2) = /'"(O) III ^^ {u + 3 M2) (w - U2f. (5.16) 

In view of Remark 2.6 the threshold function A, associated with /,, b*, c*, and 

By Theorem l5.6[ for the value a* := A;(m* ) there exists also a traveling wave tra- 
jectory connecting u* to u*^ :- U2, which we denote by v* - v*{u). By definition, 
in the phase plane it satisfies 

v'^{u) + a'^v'{u) = G:{u), (5.18) 

au Cj(m) 

with 

G'M = (mu) - /.(«*) - /.'("2) {u - ig) ^. 

We consider also the traveling wave trajectory u h^ v = v{u) connecting Mq 
to U2 which is associated with the data /, b, c\, and C2 and the threshold value 
a := A'^{iio). We will now establish lower and upper bounds on A^(mo); see (|5.241 
and (|5.251 below. 

Case 1 : First of all, in the easy case that A^(uq) (1 - e) < A;(u* ), we immediately 
obtain by l|5.17l l and then ||5.14|| 

A\uo) < {l + 2e)AHul) = (l + le) ^/^k'-^^^^u; 
.(1.2.)V3^M„o(l..) 
.(l.OVsT^Muo, 
which is the desired upper bound for the threshold function. 

Case 2 : Now, assume that A^(uo) (1 - £) > ^.(Wg) and let us derive a similar 
inequality on A''(mo)- Since G'(m2) = G',{u2) = 0, G"(m2) = G:'(m2) = 0, and 



dv dv* 

du ' du 



v{u2) = v{iQ = 0, — (M2) < 0, -j::{u2) < 0, 
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it follows from the equation 

dv b{u) G;,(m;mo,A) 

— (m + a —— = 

du ci(u) v(u) 

by letting u ^ U2 that 

This tells us that in a neighborhood of the point M2 the curve i? is locally below 
the curve v* . 

Suppose that the two trajectories meet for the "first time" at some point 

U3 e (m2, uq], so 

v{u3) = V (U3) with ^("3) > "^("3)- 
From the equations H5.5|l satisfied by J3 = v{u) and u* = v*{u), we deduce 



^ 1^(1(3)^ + a I v{u) ^^ dw = G(m3) - G(m2), 
2 J„, ci(m) 



and 



^ z;*(m3)' + a' {' v'{u) -^ du = ^(Ma) - G.(m2), 

respectively. Subtracting these two equations and using (|5.13|l and | |5.16|| , we 
obtain 

= G(m3) - G(m2) - (g.(m3) - G,{U2)) ^^-^^^ 

> (0(Mo) - C £) (M3 + 3 M2) ("3 - U2f. 

But, by assumption the curve v is locally below the curve v* so that the left-hand 
side of l|5.191 is negative, while its right-hand side of l|5.191 is positive if one 
chooses Mo sufficiently small. We conclude that the two trajectories intersect 
only at U2, which implies that u* < Uq and thus 

\v{u)\du> \v{u)\du. (5.20) 

Ju2 ^112 

On the other hand we have by i5.14|) 

A\uo)^Al-£) f"\v{u)\du<A\uo) f°^\v{u)\du ^^^^^ 

ci(0) J„, J„, ci(m) (5.21) 

= G(M2) - G(mo). 
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Now, in view of the property (i) in (|5.14l l we have 

Mo /, „ s , , 1 + 3£ , , Mn /^ . , , 3 + £ 

13 "2 + Mol < y (1 + 3 £) < \U2\ Y^' l"2 - "Ol < y (3 + £) < \U2\ Y^. 

Based on these inequalities we deduce from | |5.13| | that 

G(W2) - G(Mo) < /"'(O) ^ ^ (1 + Ce). (5.22) 

Ci((J) 

Concerning the second curve, v* = v*{ii), we have 

AS(Mo) ^ r " l^*(")l du = G,(M2) - G,(M*) 



"2 

|4 



.,,«,||^a.., 



(5.23) 



by using l|5.16t . 

Finally, combining l |5.20|l - l|5.23|l we conclude that for every £ and for all 
sufficiently small uq: 

A^(mo)<(1 + C£)aS(Wo) 

which is the desired upper bound. Exactly the same analysis as before but 
based on the cubic function /,(«) = ku^ with k = {1 - £)/"'(0) (exchanging the 
role played by /, and /, however) we can also derive the following inequality 

A\uo) > VST^ ^i^ip MO (1 - C £). (5.25) 

The proof of Theorem 15.71 is thus completed since £ is arbitrary in ( |5.24t and 
ll525ll . D 

6 Traveling waves corresponding to a given diffusion- 
dispersion ratio 



Fixing the parameter a, we can now complete the proof of Theorem 14.31 by 
identifying the set of right-hand state attainable from uq by classical trajectories. 
We rely here mainly on Theorem 15.11 (existence of the nonclassical trajectories) 
and Theorem l5.5l (critical function). 

Given mq > and a > 0, a classical traveling wave must connect m_ = uq to 
u+ = Ml for some shock speed A e {A\uo),f'{uQ)). According to Theorem 15.51 
to each pair of states {uq, u-i) we can associate the critical ratio A{uq, u-i). Equiv- 
alently, to each left-hand state Mq arid each speed A, we can associate a critical 
value B{A, Uq) = A{uq, U2). The mapping 

A 1-^ B(A, Mo) 
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is defined and decreasing from the interval A''(mo), Ao(mo) onto p,A''(Mo) ■ It 
admits an inverse 

a ^ Aa(Mo), 

defined from the interval 0, A^ (mq) onto A'' (mq)/ Ao(wo) ■ By construction, given 
any a e (0,A^(mo)) there exists a nonclassical traveling trajectory (associated 
with the shock speed A„(mo)) leaving from Mq and solving the equation with the 
prescribed value a. 

It is natural to extend the definition of the function A^ (mq) to arbitrary values 
a by setting 

A„(mo) = A''(mo), a>A\uo). 

The nonclassical traveling waves are considered here when a is a fixed param- 
eter. So, we define the kinetic function for nonclassical shocks, 

(mo,«)i-^ (piiuo) = "2, 
where 112 denotes the right-hand state of the nonclassical trajectory, so that 

Mo -M2 

Note that <Pa{uo) makes sense for all uq > but a < A^(uo). 
Theorem 6.1. For all Uq> and a > and for every speed satisfying 

Aaiuo) < A <f'(uo), 

there exists a unique traveling wave connecting w_ = uq to u+ = ui. Moreover, for 
a > A^{uq) there exists a traveling wave connecting U- = uq to u+ = uifor all 

Ae[A^(wo),/'(wo)]. 

Proof. We first treat the case a < /^{uq) and A e [Aa{uQ),f'{uQ)\. Consider the 
curve u 1-^ V-(u; A, a) defined on [mi, Mq] that was introduced earlier in the proof 
of Theorem l5.ll We have either C-(mi; A, a) = and the proof is completed, or 
else C-(mi; A, a) < 0. In the latter case, the function V- is a solution of HS.SI I that 
extends further on the left-hand side of U- in the phase plane. On the other 
hand, this curve cannot cross the nonclassical trajectory u 1-^ v{u) connecting 
U- = Uq to u+ = (pJ',.(mo)- Indeed, by Lemma l46l we have 

Ji(uo;A,a) < Jl{uo;Aa{uo),a). 

If the two curves would cross, there would exist u* e ((pj^(uo)/ Mi) such that 

/ .^ / ,x 1 dv , ,^ dv- , ,, 

v(u) = v-{u) and — (w ) < ^— (« )• 
du du 
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By comparing the equations (|5.5l l satisfied by these two trajectories we get 



This leads to a contradiction since the right-hand side of (|6.2|l is positive while 
the left-hand side is negative. We conclude that the function V- must cross the 
M-axis at some point ut. with 112 < (Pa(wo) < "3 < mi- The curve u 1-^ i'-(m. A, a) on 
the interval [M3, Mq] corresponds to a solution y 1-^ m(i/) in some interval (-00, 1/3] 
with Uy{y^) = and 

^ _ g{u{y3),A)-g{uo,A) _ g;,(»3;»o,A) , , „. 

which is positive by Theorem l4.7l ThusMyy(i/3) > and necessarily m(i/) > M3for 
y > yo,. Indeed, assume that there exists 1/4 > 1/3, such that 1/(1/4) = "(ys) = "3- 
Then, multiplying | |4.25l l by V-IC2 and integrating over [1/3, 1/4], we obtain 

1 . p &(m) . 

7:i^_(i/4) + a —r\ — ;— - 1^_ dy = G(m3; uo. A) - G(m3; Mo, A) = 0. 

-^ Jy3 ClWC2\U) 

This would means that u{y) = 113 for all y, which is excluded since U- = U\ . 

Now, since u < uq we see that u is bounded. Finally, by integration over the 
interval (-00, y] we obtain 

1 , ry b(u) , 

which implies that v is bounded and that the function u is defined on the whole 
real line R. When y ^ +00 the trajectory (m, c) converges to a critical point 
which can only be (mi, 0). 

Consider now the case a > A^{uq). The proof is essentially same as the one 
given above. However, we replace the nonclassical trajectory with the curve 
u h^ v+{u) defined on the interval [u2/Mi]- For each A fixed in (A^(uo)//'(wo)) 
(since a > A^{uo)) and thanks to Remark 4.4, the function, W = V+-V- (defined 
in the proof of Theorem 15. 1[ with V-{u;A,a) and v+{u;\,a) and extended to 
A e (/'(m2),/'(wo))) satisfies W(a) < 0. On the left-hand side of mj, with the 
same argument as in the first part above, we can prove that the extension of V- 
does not intersect v+ and must converge to (mi,0). Finally, the case A = A^{uq) 
is reached by continuity. This completes the proof of Theorem 16. II D 

Theorem 6.2. If A^{uq) < A < Aa(Mo) there is no traveling wave connecting u_ = Uq 
to u+ = U\. 

Proof. Assume that there exists a traveling wave connecting mq to mi. As in 
Lemma [5.21 we prove easily that such a curve must approach (mq/O) from the 
quadrant Q\ and coincide with the function V- on the interval [mi,Mo]- On the 
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other hand, as in the proof of Theorein l6.11 we see that this curve does not cross 
the nonclassical trajectories. On the other hand, we have 

fr(Mo; A, a) > Jl{uo; Aa(Mo), «), 
thus, the classical curve remains "under" the nonclassical one. So we have 

^-(<Pa(Wo)) < v((pi{UQ)), 

where u i-^ (m, v{u)) denotes the nonclassical trajectory. Assume now that the 
curve {u,V-{u)) meets the !<-axis for the first time at some point (1^3,0) with 
M3 < (Pa(Mo) < "2- The previous curve defined on [iis, mq] corresponds to a 
solution y 1-^ u{y) defined on some interval (-00,1/3] with Uy{y^) = and 
Myyd/s) ^ 0. Thus ^1/(1/3) Satisfies (|6.3t and is negative (Lemma l5.3t . This implies 
that Uyy{y^) < which is a contradiction. Finally, the trajectory remains under 
the u-axis for u < 1/2, and cannot converge to any critical point. D 



According to Theorem 16.11 the kinetic function can now be extended to all 
values of a by setting 

cpi{uo) = cp\uo), a>A\uo). (6.4) 

Finally we have: 

Theorem 6.3. (Monotonicity of the kinetic function.) For each a> Othe mapping 
Uq i-> <Pa("o) is decreasing. 

Proof. Fix Mo > 0, a > 0, A = Aa(Mo) and 112 = <Pa{uo)- First suppose that 
a > A\uo). Then, for all w* > uq, since cp^ is known to be strictly monotone, it is 
clear that 

<pi{ul) < <p\ul) < <p\uo) = <pi{uo). 

Suppose now that a < A^{uo). Then, for m* > uq in a neighborhood of uq, the 
speed A* = — °,_^ satisfies A* e (A''(Mg), Ao(Uq)j. Then, there exists a nonclas- 
sical traveling wave cormecting u* to 112 for some a* = A{ii'Q, 112). The second 
statement in Theorem 15.51 gives a* > a. Since the function A^ is decreasing 
(by the first statement in Theorem 15.51 we have Aa,(Mg) < Art(Mg) and thus 
(PaiUg) < U2 = (Pai^o) and the proof of Theorem l6.3l is completed. n 

Proof of Theorem [Ol Section 5 provides us with the existence of nonclassical 
trajectories, while Theorems l6.1l and l6.2l are concerned with classical trajectories. 
These results prove that the shock set is given by l|4.14|l . By standard theorems 
on solutions of ordinary differential equations the kinetic function is smooth 
in the region \a < A'^^{uo)\ while it coincides with the (smooth) function cp^ 
in the region \a > A\uo)\. Additionally, by construction the kinetic function 

is continuous along a = A^(uo). This proves that cp^ is Lipschitz continuous 
on each compact interval. On the other hand, the monotonicity of the kinetic 
function is provided by Theorem 1 6. 3 1 The asymptotic behavior was the subject 
of Theorem 15. 71 □ 
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Remark 6.4. To a large extend the techniques presented in this paper extend to systems 
of equations, in particidar to a classical model of elastodynamics and phase transitions. 
The corresponding traveling wave solutions {v, zv) = {v{y), w{y)) must solve 

-S Vy - 1L{w, Wy, Wyy) = {^p{w) V y^ , 
S Wy +Vy^O, 

where s denotes the speed of the traveling wave, E is the total stress function, and p{w) 
is the viscosity coefficient. When E is given, some integration with respect to y we 
arrive at 

\'{w) T j \ 

-s {v - V-) - a{w) + a{w-) - p{w) Vy = — — Wy - [A{w) Wyj , 

S {zv - W-) + V-V- = 0, 

where {v-, zv-) denotes the upper left-hand limit and A(w) the capillarity coefficient. 
Using the second equation above we can eliminate the unknozvn v{y), namely 

A{w)^l^ i\{w)^^^ Wy\ + p{w)Vy = s^ {zv -ZV-) - a{zv) + a{w-), (6.5) 

which has precisely the structure of the equation | |4.4|I studied in the present paper, so 
that most of our results extend to the equation | |6.5| I. 
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